Connecting Lattice QCD with Chiral Perturbation Theory at Strong Coupling 
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We study the difficulties associated with detecting chiral singularities predicted by chiral pertur- 
bation theory (ChPT) in lattice QCD. We focus on the physics of the remnant 0(2) chiral symmetry 
of staggered fermions in the strong coupling limit using the recently discovered directed path algo- 
rithm. Since it is easier to look for power-like singularities as compared to logarithmic ones, our 
calculations are performed at a fixed finite temperature in the chirally broken phase. We show that 
the behavior of the chiral condensate, the pion mass and the pion decay constant, for small masses, 
are all consistent with the predictions of ChPT. However, the values of the quark masses that we 
need to demonstrate this are much smaller than those being used in dynamical QCD simulations. 
We also need to use higher order terms in the chiral expansion to fit our data. 



INTRODUCTION 

With the advent of fast computers it is now possible 
to calculate many hadronic quantities from first princi- 
ples using lattice QCD. However, today these calcula- 
tions contain systematic errors due to finite lattice spac- 
ings, finite volumes and quenching. Although in principle 
all these errors can be controlled, the only clean way to 
reduce quenching errors is to perform unquenched calcu- 
lations. Due to algorithmic difficulties today most un- 
quenched calculations are performed using large quark 
masses and the results are then extrapolated to realistic 
values using (quenched and partially quenched) ChPT. 
Unfortunately, recent attempts to connect lattice QCD 
with the usual one-loop ChPT predictions have failed to 
give clear answers 0, S S • It is now believed that lat- 
tice artifacts should be taken into consideration in ChPT 
0- 0- 0- 01 ■ It has been suggested in 0,^3 the lattice data 
is described better by the resulting more elaborate fit- 
ting functions. There are also other interesting attempts 
to extract useful information using finite size effects in 

ChPT [ni. 

Given the difficulties associated with understanding 
chiral singularities in a realistic calculation of QCD, in 
this paper we explore the subject in strong coupling lat- 
tice QCD with staggered fermions. We use a very effi- 
cient algorithm discovered recently to solve this model in 
the chiral limit Although the strong coupling limit 
suffers from severe lattice artifacts, when the quarks are 
massless lattice QCD with staggered fermions has an ex- 
act 0(2) chiral symmetry which is broken spontaneously. 
Thus, our model contains some of the remnant physics of 
chiral singularities expected in QCD. In particular, there 
are light pions and it would be useful to understand the 
range of quark masses where the singularities predicted 
by conventional ChPT can be seen. 

Instead of focusing on chiral singularities that are loga- 
rithmic, in this work we focus on power-like singularities 
that arise at finite temperatures. Recently, it was shown 
with high precision that our model undergoes a second 
order chiral phase transition at a critical temperature 



T r . This transition belongs to the 0(2) universality class 
Thus, at a fixed temperature below T c , within the 
the 0{2) scaling window, the long distance physics of our 
model is described by a three-dimensional 0(2) field the- 
ory in its broken phase. At this temperature there is a 
range of quark masses where the light pions are describ- 
able by the conventional continuum ChPT. The effective 
action can be written in terms of 5, an 0(2) vector field 
with the constraint 5-5=1. At the lowest order this is 
given by 0] 
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where h is the magnetic field and \h\ is identified with the 
quark mass. The two low energy constants F 2 and £ that 
appear at this order are the pion decay constant and the 
chiral condensate respectively. We note that since we are 
discussing a three-dimensional effective theory, F 2 has 
the dimensions of inverse length. Thus, we can define 
a correlation length £ = 1/F 2 . For massless quarks £ 
diverges as (T c — T)~ v at T c . In order to connect our 
model with ChPT and avoid lattice artifacts we need 
to choose the temperature and quark masses such that 
1 << £ << l/M^ (lengths are measured in lattice units) 
where M T is the pion mass. In this region we expect 
observables such as the chiral condensate, the pion mass 
and the pion decay constant satisfy the expansion 



(O) = z Q + zi^/m, + z 2 m + z 3 m^/m + 



(2) 



where the y/m behavior is the power-like singularity aris- 
ing due to the infrared pion physics [lfij . The main ob- 
jective of this paper is to detect these singularities in 
the context of strong coupling lattice QCD with stag- 
gered fermions. Such power-like singularities have been 
observed in spin models jig , but as far as we know have 
not been studied with precision in lattice QCD calcula- 
tions. 
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THE MODEL AND OBSERVABLES 

The partition function of the model we study in this 
article is given by 

Z(T, m) = J [dU] exp (S[U, ip, $]) , (3) 

where [dU] is the Haar measure over U(3) matrices and 
[dip dip] specify Grassmann integration. At strong cou- 
pling, the Euclidean space action S[U, ip, ip] is given by 



2 

X,fl 



m ip x ip x 



(4) 

where x refers to the lattice site on a periodic four- 
dimensional hyper-cubic lattice of size L along the three 
spatial directions and size L t along the euclidean time 
direction. The index fj, — 1,2,3,4 refers to the four 
space-time directions, U x ,n € U(3) is the usual links 
matrix representing the gauge fields, and ip x ,tp x are t ne 
three-component staggered quark fields. The gauge fields 
satisfy periodic boundary conditions while the quark 
fields satisfy either periodic or anti-periodic boundary 
conditions. The factors T) x „ are the well-known stag- 
gered fermion phase factors. We choose them to be 
rjgn = = 1,2,3 (spatial directions) and 77^4 = T 
(temporal direction), where the real parameter T acts 
like a temperature. By working on asymmetric lattices 
with Lt < < L at fixed Lt and varying T continuously one 
can study finite temperature phase transitions in strong 
coupling QCD |l7j. We use U(3) gauge fields instead of 
SU(3) in order to avoid inefficiencies in the algorithm due 
to the existence of baryonic loops. This distinction is not 
important for our study since the baryons are expected 
to have a mass close to the cutoff. 

The partition function given in Eq. J2J can be rewritten 
in a monomer-dimer representation as discussed in detail 
in 0,0. Every configuration in the new representation 
is described by monomer variables n x = 0, 1, 2, 3 associ- 
ated to the sites, and dimer variables b Xtll — 0,1,2,3 
associated to the bonds connecting neighboring sites x 
and x + ft, along with the constraint that at each site, 

n x + X^x,^ + K-hn] = 3. 

In order to study chiral physics we focus on the follow- 
ing observables: 



(i) The chiral condensate 



(5) 



(ii) the chiral susceptibility 
1 1 d 2 



X 



L 3 Z dm 2 



Z(T,m), 



(6) 



(iii) the helicity modulus 

Y = f ({E J ^ 2 + E J ^ 2 + E J M] a }). (7) 

\ X X X f 

where Jx,^ — <^x{^>x,u ~ N/8), with a x — 1 on even 
sites and ctj = -1 on odd sites and 

(iv) the pion mass, obtained using the exponential de- 
cay of the correlation function along one of the spa- 
tial directions fi — 1,2,3: 

lim V(a x n x n )=Ce- M -l^l, (8) 

IzJ^oo — 

x ± 

where x± refers to components of the coordinate x per- 
pendicular to the n direction. When m = the cur- 
rent Jx.fi is the conserved current associated with the 
0(2) chiral symmetry. As discussed in [14|, one can 
define the pion decay constant at a quark mass m to 
be = \imL_tooY. For m = we then obtain 
F = lim m ^o -Pmi the pion decay constant introduced in 
Eq.JI). We can also define S m — limi_ >oo (0) the infinite 
volume chiral condensate. Again E = lim m _ >o S m . 

RESULTS 

We have done extensive computations at L t = 4. It 
has been shown with high precision in 13] that this 
model undergoes a chiral phase transition which belongs 
to the three-dimensional 0(2) universality class at a criti- 
cal temperature T c — 7.47739(3). FigureQlshows the plot 
of F 2 as a function of T. Based on universality we expect 

S~ =F 2 = C(T C ~ T) v for T <T C with v = 0.67155(27) 
If we fit our results to this form in the range 
7.05 < T < 7.42, with T c and v fixed to the expected 
values, we find C = 0.2217(1) with x 2 /DOF= 1.03. In- 
cluding T = 7.0 in the fit makes the % 2 /DOF jump to 
2.0. At a temperature very close to T c , we expect £ to 
be extremely large and it would be difficult to satisfy 
l/Mjr >> £ with our limited computing resources. On 
the other hand in order to avoid lattice artifacts it is im- 
portant not to have £ ~ 1. Using the above analysis we 
estimate that T = 7.0 is at the edge of the 0(2) scaling 
window and hence we choose to fix T at this value and 
study the long distance physics near the chiral limit. Wc 
vary the spatial lattice size from L = 8 to L = 96 for a 
variety of masses < m < 0.025. At m = we also have 
data at L = 120, 144. Below we discuss our main results. 

The Pion Decay Constant 

Our results show that for L > 32 the finite size effects 
on Y are smaller than the statistical errors. Hence it is 
relatively easy to compute F^ by extrapolating results 
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for Y at different volumes. In Table[I]we give the results 
for some values of the quark mass. The first four terms 



m 




m 




0.00000 


0.3675(1) 


0.0050 


0.4063(1) 


0.00025 


0.3703(1) 


0.0075 


0.4188(1) 


0.00050 


0.37300(15) 


0.0100 


0.42890(15) 


0.00100 


0.37804(15) 


0.0200 


0.4583(2) 


0.0025 


0.3903(2) 


0.0250 


0.4693(2) 



TABLE I: F m as a function of m. 

in the chiral expansion of F m in three-dimensional 0(N) 
ChPT are given by 14] 



F m = F[l + a\\pm + a,2m + a^m\/rn] 1 



(9) 



where a\ oc (N — 2). Since in our case TV = 2 we expect 
ai = 0. In order to check if our results fit the predic- 
tions of ChPT we fit our data to Eq.@ in the range 
< m < 0.00625. We find F = 0.3674(1), ai = 0.01(2), 
a 2 = 34.0(7) and a 3 = -182(6) with x 2 /DOF= I- 2 - Our 
data and the fit are shown in Fig. [21 The prediction 
of ChPT that a\ = is in excellent agreement with our 
results. Fixing a x = in the fit yields F = 0.36747(6), 
a 2 = 34.3(3) and a 3 = -185(3), while the x 2 /DOF re- 
mains essentially unchanged. Interestingly if we fit the 
data in the range < m < 0.00135 with a\ = a 3 = 
fixed, we find that F = 0.36764(7) and a 2 = 28.3(3) with 
X 2 /DOF= 0.9. This shows that there are systematic er- 
rors in evaluating the fitting parameters due to contami- 
nation from higher order terms. In the inset of Fig. [2]we 
focus on the extremely small mass region and show that 
ai = is even clear to the eye. This is one of the main 
results of our paper. We suggest that ai = is a useful 
signature of the 0{2) universality and could be used in 
future studies of lattice QCD with staggered fermions. 



The Chiral Condensate 

Chiral perturbation theory predicts that the first four 
terms in the chiral expansion of E m are given by 



S[l + bi\/m + b 2 m + b^my/m] 



(10) 



As shown in [14|. b\ oc (N — 1), implying that b\ ^ in 
our case. In order to test how well Eq. I|10|l describes 
our model, we compute S m at small masses. 

At m = we compute £ by using the finite size scaling 
formula for \ given by (l4l | 



X : 



-j: 2 l 3 

N 



1 + 0.226(^-1)^ + ^ 



(11) 



where a depends on the higher order low energy con- 
stants of the chiral Lagrangian. When F is fixed to 0.3675 
obtained earlier, our results for 16 < L < 144 fit very well 
to this formula. We find £ = 2.2648(10), a = 4.6(3) with 
X 2 /DOF= 0.87. 



m 




m 




0.00000 


2.2648(10) 


0.0050 


2.6340(08) 


0.00025 


2.2340(10) 


0.0075 


2.7278(15) 


0.00050 


2.3560(10) 


0.0100 


2.8025(15) 


0.00100 


2.4040(05) 


0.0200 


3.0199(11) 


0.0025 


2.5103(08) 


0.0250 


3.0986(08) 



TABLE II: Values of E m at selected values of m. 

For m ^ one can measure £ m by extrapolating the 
results for ((f) at different volumes. However, due to crit- 
ical slowing down it is difficult to measure the condensate 
accurately at a fixed small mass on large volumes. With 
our computing resources we find that this procedure gives 
reliable answers for £ TO only up to m ~ 0.001. On the 
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FIG. 3: Plot of E m vs. m. The solid line represents the FI G- 4: Plot of M w vs. m. The solid line represents the 
function E m = 2.2642[1 + 1.36Vro + 23m - 135m Vm]. function M n = 4.095Vm[l + 0.68Vm - 25.8m + 149mVrn]. 



The Pion Mass 



other hand when to ^ 0, in the large volume limit we 
expect x ~ ^m-^ 3 U P to an additive constant plus ex- 
ponentially small corrections. We find that the signal 
for S m extracted from fitting the data to this form is 
much cleaner at small masses than the signal obtained 
from the direct measurement. Further, the values of 
E m obtained using this procedure agree well with the di- 
rect measurement at larger quark masses. In table ITT1 we 
give the values of £ m obtained for selected values of the 
quark masses. Fitting our results to Ea. flop in the region 
< m < 0.00625 gives S = 2.2642(10), h = 1.36(4), 
b 2 = 23(1), 6 3 = -135(10) with a x 2 /DOF= 1.1. Our 
results are shown in Fig [31 We note that we cannot find 
a mass range within our results in which we can find 
a good fit when we fix 62 = &3 = 0. However, in the 
range < to < 0.00175 we can set 63 = to obtain 
S = 2.2643(10), 61 = 1.44(3) and b 2 = 15.2(7) with a 
X 2 /DOF= 1. Note that b 2 changes by about 30% when 
this different fitting procedure is used, while 61 is more 
stable. 



The first four terms in the chiral behavior of the pion 
mass are predicted to be of the form 



1 + c\ \pm + C2TO + c^ra^/m 



(12) 



Further, the chiral Ward identities imply that c\ — bi/2 — 
ai. Our estimates for the pion masses for a selected range 
of m are shown in Table UTTl Fitting our data to Eq. i|12|) 
after fixing £ = 2.2648 and F = 0.3675 obtained above, 
we find d = 0.55(16), c 2 = -21(6) and c 3 = 110(50) 
with x 2 /DOF= 0.5. We see that our data is consistent 
with the relation c± = bi/2 — ai although the error in c\ is 
large. Fixing c\ = bi/2 — 0.68 obtained from fitting the 
chiral condensate, yields c 2 = —25.8(8) and C3 = 149(14) 
without changing the quality of the fit. This latter fit 
along with our results for M n are shown in Fig. 



DISCUSSION 



m 




m 




0.00000 


0.000 


0.0050 


0.2812(8) 


0.00025 


0.0650(3) 


0.0075 


0.3399(5) 


0.00050 


0.0920(4) 


0.0100 


0.3871(4) 


0.00100 


0.1295(4) 


0.0200 


0.5285(10) 


0.0025 


0.2022(3) 


0.0250 


0.5837(8) 



TABLE III: Values of M n at selected values of m. 



It is expected that the natural expansion parameter for 
ChPT in three dimensions is x = A/ 7r /(47rF 2 ). Let us find 
the values of x where the chiral expansion up to a certain 
power of x is sufficient to describe the data reasonably. 
Of course this question is model-dependent, since some 
models may have larger contributions at higher orders 
compared to others. Here we ask this question in the 
context of the model studied in this paper. 

In order to understand the dependence of and E m 
on x we plot in Fig. these two observables as func- 
tions of x. The solid lines are the best fits to the form 
zq + z\x + z 2 x 2 + z^x 3 , while the dot-dashed lines are 
the best fits to a smaller range in x with Z3 set to 0. As 
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FIG. 5: Plot of and E m vs. x = M„/(4tyF 2 ). The solid 
lines and dot-dashed lines are fits to two different orders of 
chiral expansion discussed in the text. 

the graph indicates the error in E m due to the absence 
of the Z3 term is different from the error in at a given 
value of x. In order to determine within say 5% we 
need the Z3 term even at x ~ 0.15. This shows that our 
model contains important higher order terms. An inter- 
esting question which we cannot answer at this point is 
whether this property is generic or not. In any case we 
have shown that the connection of lattice QCD data to 
ChPT is indeed possible as long as one does not assume 
that higher order terms are negligible. For the lowest or- 
der terms to dominate it may be necessary to go to much 
smaller quark mases than has been possible until now. In 
fact it is much easier to connect our results with ChPT at 
?7i = 0, which is called the e-regime of ChPT [l^. Thus, 
finding an algorithm to work directly at m — is a useful 
goal to strive for in the future. Finally we hope that our 
results will motivate further work to uncover the power- 
law chiral singularities at finite temperatures in Lattice 
QCD at weaker couplings. This should be easier than 
looking for logarithmic singularities at zero temperature. 
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